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Introduction In a 2001 article, Fetecau and Zierep f\\ considered Stokes' first problem for a second grade 
(SG) fluid, unaware that it had already been solved (correctly) by Puri |2| in 1984. These authors used the 
Fourier sine transform to obtain, what they believed to be, the exact solution of the resulting initial-boundary- 
^ . value problem (IB VP). Unfortunately, due to their incorrect computation of the distributional derivative of the 

Heaviside function, an omission of a critical term in the subsidiary equation occurred, an elementary mistake 
(gracefully explained in |3|) apparently not uncommon in the literature. Therefore, the solution given in [I] 
^ ■ is, generally speaking, incorrect. 

(^_| \ The aim of the present comment is to point out and correct this error. Hence, employing the same notation 

^ • convention used in llj, we begin by restating the IB VP considered therein: 



(V + ad,)djv(y, t) = d,v{y, t), {y, t) e (0, cx,) x (0, cx,); (la) 
v(j, 0) = 0, y>0; v(0, f) = V{t), v(y, t), d,-v{y, f) ^ as y ^ co, f > 0. (lb) 



Here, a and v are positive constants. Note that V{t) as written above and in 01 is an ambiguous notation 
because the function is not arbitrary. For this IB VP to be Stokes' first problem, the plate must move suddenly 
' at f = 0^. Therefore, V'(f) must be such that V(t) - V(t)6(t), where 0{-) denotes the Heaviside unit step function 

and V{t) is some smooth function. (Stokes himself made this distinction clear [4 pp. 101].) Consequently, the 
correct expression in lieu of V'{t) in Uj Eq. (2.3)] is V'{t)6{t) + V(t)6(t), where 6(-) denotes the Dirac delta 
function. In the case of V(t) ~ as ? ^ 0, for some integer « > 1, the singularity is ameliorated and by pure 
i luck the proper answer may be obtained despite the erroneous derivation. 

I Corrected solution The simplest case is that of «m/)MZi;ve plate motion with constant velocity: v(0, f) - Vo9{t). 

. Thus, V'{t) in I becomes Vo6{t) and it follows that the solution of the subsidiary equation (i.e., H] Eq. (2.3)]) 

' is, in this case, given by 



v,(^, s) = 2Vo \ — — — , + — — — ^ , \ , (2) 



H I which we obtained using the Laplace transform. Here, a bar over a quantity denotes its image in the Laplace 

^ transform domain and s denotes the Laplace transform parameter. Note that we employ the Fourier sine 

transform convention from IJl. 
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Using partial fractions, a standard table of Laplace inverses, and the definition of the inverse sine trans- 
form, Eq. (|2]l is easily inverted to yield the time-domain solution 



nJa ^ ^\l+ae) ^ n Jo l+< ^\i+aej\ 



(3) 



Here, we observe that the second integral in Eq. (O, which represents the contribution from the second quo- 
tient in Eq. (|2|, is «of present in fT, Eq. (2.6)]; consequently, the latter is incorrect. This error appears to be due 
to the authors of |1 1 taking V'{t) as zero, rather than its correct value VoSit). Equation Q above is, therefore, 
the correct exact solution, based on the Fourier sine and Laplace transforms, to Stokes' first problem for a SG 
fluid. 

The same mistake is committed when going from Eq. (3.3) to Eq. (3.4) in [1], however, we leave it as an 
exercise to the reader to find the correct solution in this case. 

An additional error transpires in the integration of the subsidiary equation \\, Eq. (2.3)]. Though it is 
claimed that IT] Eq. (2.4)] is its solution, it is clear that, by the integrating factor method, which appears to 
be what was attempted, such a solution cannot be obtained. The same mistake occurs in the derivation of 11] 
Eq. (3.3)]. We find it unnecessary to supply corrected expressions for these intermediate formulae because 
they already suffer the much graver error, committed in the taking the derivative of the Heaviside function, 
that we corrected above. 

A simple finite-difference scheme To provide an independent check on the corrected transform solution given 
in Eq. ([3]l, we also solve the IB VP ^ numerically. We define the (uniform) spatial- and temporal-step sizes 
Ay := L/(M-l) and At := tf/{K-l), where M > 2 and A' > 2 are integers and now (j, f) e (0, L)x(0, f/]. Then, 
letting d" ^ v(jj, t") be the approximation to the exact solution on the grid, where yj :- jAy (0 < j < M - 1) 
and t" :- nAt (0 < n < A" - 1), we may discretize Eq. (flal) as follows: 



6,^v'} - aS,^6,^6y^v'} = yl [dy^Sy^v"/' + 5,+5v-C ") ■ (4) 

Here, St+ is the forward temporal difference operator and 6y+ and Sy- are, respectively, the forward and back- 
ward spatial difference operators |6 §3.3]. It is a straightforward, though lengthy, calculation (see, e.g., ||6l) 
to show that this implicit, two-level Crank-Nicolson-type discretization is unconditionally stable and has 
truncation error 0[{At)^ + (Ay)^]. The boundary conditions are implemented as 

"-{vo, Ill'^K-ir (5) 

For appropriately chosen L » 1, the front does not reach the y - L boundary for any / e (0, //], and so this is 
the "numerical infinity." 

For the computations shown below, we took M - 5000, K - 1000 and L = 20 to obtain a highly-accurate 
solution that we can compare to Eq. (O. Matlab's built-in Gaussian elimination algorithm was used to invert 
the symmetric tridiagonal matrix resulting from this discretization. 

Illustrated example Here we compare the correct solution given in Eq. (|3]l and the numerical solution by the 
difference scheme in Eq. (|4]l to the incorrect solution given in [T, Eq. (2.6)]. Clearly, Eq. ([3]l and the numerical 
solution are in excellent agreement. Meanwhile, the solution from IT] Eq. (2.6)] does not even satisfy the 
boundary condition lim,,^o+ v^y, f) = Vq (? > 0). 



Conclusion It is regrettable that this error has now spread throughout the literature on non-Newtonian fluid 
mechanics, "infecting" a number of journals over the last nine years. The same error also appears in work on 
Stokes' first problem for relaxing fluids |7 ]. This means that one must be very careful when applying repeated 
integral transforms to problems with mixed derivatives and discontinuous initial and/or boundary conditions. 

It is curious to note that if the contrived source term 2VQa5{t)5' iy) were added to the left-hand side of 
Eq. ( [Tal l, the error in computing V'(f) would be "cancelled" and IT] Eq. (2.6)] would be the solution to this 
new problem. Similarly, if the boundary condition in Eq. dlbb were changed to v(0, f) = (1 - &^^'^")VQ0{t), 
then [1 Eq. (2.6)] coincides with the well-know solution to this problem from the theory of unsteady flow 
in fissured rocks [8 Eq. (5.9)]. This is the kind of lucky occurrence we alluded to in the Introduction. It is 
noteworthy that [8] is a seminal paper with several hundred citations, yet it somehow eluded the authors of 

m. 



3 



v/Vo 





Fig. 1 V vs. >> for (left) v = 0.l,a = 0.01, t = 0.01 and (right) v = 0.1, o- = 0.1, f = 0.5. Dashed: incorrect solution (1 Eq. (2.6)]; 
solid: correct solution [Eq. above]; dots: numerical solution [using the scheme in Eq. Q]. 



In addition, it is of interest to note that if we set 



2 f 

u - v/Vo, Y -y, T - vt, I - a, rj - 



then Eq. (O can be re-expressed as 

m(f, t) = e{T) 
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(6) 



(7) 



The latter is equivalent to the b - Q special case of fT, Eq. (10)], wherein the correct porous medium version 
of IB VP ^ is solved, a result which can also be obtained using only the Laplace transform with respect to 
time 119 2 1. This means that it is possible to prevent the error corrected herein by avoiding the use of multiple 
integral transforms simultaneously. 

Finally, the assumption that a is positive in Eq. ( [Tal l is an important one because a < Q renders the 
rest state of the SG fluid unstable [10] and the problem ill-posed. Yet, certain theoretical arguments and 
experimental measurements involving polymer solutions suggest that a may be negative llll . All this serves 
to point out the approximate nature of the SG fluid model |12| . In fact, the sign of a was not addressed in ||T||, 
and the correction of the mathematical solution we have presented should not be construed as resolving any 
fundamental difficulties with the mechanics of SG fluids. 
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